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1. INTRODUCTION 


CuRISTIAAN HuyGeNns! noticed that one of his specimens of iceland spar 
contained an inner layer exhibiting iridescence and remarked that this was 
parallel to one edge of the rhombohedral cleavage and equally inclined to 
the two others. It is now recognised that such layers arise from internal 
twinning and that they may be artificially produced by subjecting a crystal 
of calcite to suitable stresses. It is also well-known that when a distant 
luminous object is viewed through a rhomb of iceland spar in which such 
twinning is present, three images of the object are seen instead of only one; 
the two outer images often exhibit vivid colour, while the central or un- 
deviated image is usually colourless. The brightness of the images, their 
angular separation, and the colours which they exhibit vary greatly with 
the direction in the crystal along which the source of light is viewed. An 
explanation in general terms of the origin of the two outer images and of 
the colour and polarisation which they display appears in Mascart’s treatise.” 
Grailich? and later Osthoff* investigated the reflection and refraction of 
light by the twinning planes in calcite. A summary of their results is given 
by Pockels.® A noteworthy observation regarding the optical behaviour 
of twinned calcite is also to be found recorded in a paper by Rayleigh.® 


Numerous specimens of transparent calcite exhibiting internal twin- 
ning were available to us in the Museum of this Institute and this induced 
us to undertake the study of their optical behaviour. It was immediately 
noticed that specimens in which twinning layers traverse the crystal simul- 
taneously in two or even three different directions exhibit gorgeous arrays 
of multi-coloured images arranged in geometric patterns when a source of 
light is viewed through them. Plates I and II accompanying this paper 
reproduce photographs obtained by us of a few of the many patterns observed 
with the material at our disposal. A detailed examination of the case of 
calcite containing only a single twinning layer also revealed many facts which 
had apparently been overlooked by the earlier investigators. We shall refer 
to these in the course of the paper, the purpose of which is to present the 
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complete picture of the optical behaviour of iridescent calcite which has 
emerged from our studies. 


2. GEOMETRIC THEORY 


When a train of light-waves traverses a layer of calcite differing from the 
material on either side of it in crystallographic orientation, the waves would 
necessarily be reflected and refracted by the layer. Elementary considerations 
indicate that the angle of incidence @; and the angle of reflection or refrac- 
tion #, would be related to the corresponding refractive indices p; and p, 
by the general formula 


Hy sin 0; = tpSin 0+. (1) 


Since the medium is birefringent, it follows that for a given angle of inci- 
dence, a wave-train of either the ordinary or the extraordinary species would, 
in general, give rise to fwo reflected and two refracted wave-trains. Other 
inferences from the formula are: (1) One of the refracted wave-trains emerg- 
ing from the layer in each case would travel in the same direction as the 
incident wave-train; (2) the angles of reflection and of refraction of the 
ordinary wave-train would necessarily be identical; and (3) the angles of 
reflection as also of refraction would, in general, be greater for the extra- 
ordinary than for the ordinary wave-trains (see Fig. | and Fig. 2 below in 
which these results are embodied). 


oO 


0 





Fic. i. Fic. 2. 
Reflection and Refraction by a Twinning Layer. 


When a distant source of light is viewed through the opposing pair of 
faces of a cleavage rhomb of iceland spar, we observe only a single image 
of the source, though the wave-trains divide into the ordinary and extra- 
ordinary species and follow different paths within the crystal. It is readily 
seen from Figs. | and 2 that the same situation would subsist even when a 
twinning layer is present and is traversed by these waves. In such a case, 
however, the ordinary wave-train which emerges as an extraordinary one 
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and the extraordinary wave-train which emerges as an ordinary one after 
passage through the layer, would both give rise to additional images, deviated 
respectively to one side and the other of the principal image. The magni- 
tude of their angular deviations may be deduced at once from the general 
formula and is the smaller, the more nearly the direction of observation 
through the crystal approaches the optic axis. It is important to remark 
that while the central or undeviated image would appear in the same direc- 
tion for all wavelengths, such would not be the case for the deviated images. 
The dispersion of the refractive indices in calcite is very different for the 
ordinary and the extraordinary waves. Hence, the two lateral images would 
appear drawn out into spectra, and it may be remarked that such dispersion 
would be large when the wave-train under consideration emerges nearly 
parallel to the twinning layer. 

The general formula indicates other results of interest. If pu; be 
the ordinary and jp, the extraordinary index, @, would be real only if 
#;< sin-'y,/y;. At the upper limit indicated, the reflected and refracted extra- 
ordinary waves would emerge parallel to the twinning layer, while for greater 
incidences they are non-existent. Likewise, when p; is the extraordinary 
and yu, the ordinary index, @,<sin“y;/u,: the limiting value is reached 
when the incident waves graze the twinning layer. In both cases, the limit- 
ing angle may be written as sin-!e/u9 where ue is the ordinary index and 
te is the extraordinary index for the direction of propagation parallel to 
the twinning layer in the plane of incidence. Thus, if light be incident on 
a twinning layer simultaneously in all directions, a boundary line may be 
drawn in the field separating the directions of incidence of the ordinary wave- 
trains which give rise to extraordinary reflections and refractions and of 
those which do not. Likewise, the optical fields traversed by the reflected 
and refracted wave-trains are, each of them, divided by a boundary line 
on the two sides of which ordinary wave-irains resulting from incident extra- 
ordinary ones respectively appear and do not appear. 


It is evident from the foregoing remarks that we cannot always expect 
to be able to observe the deviated images on both sides of the principal 
image when a distant light source is viewed through a calcite rhomb. As 
the angle of incidence on the twinning layer approaches and ultimately 
exceeds the limiting value, the extra-ordinary refraction resulting therefrom 
would emerge parallel to the twinning layer and ultimately disappear. 


3. INTENSITY CONSIDERATIONS 


A factor which greatly influences the observed phenomena is the azi- 
‘muth of the plane of incidence of the light on the twinning layer. It has 
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already been remarked that the layer is parallel to one of the edges of the 
rhombohedral cleavage and equally inclined to the two others. It follows 
that a plane through the edge first mentioned which is perpendicular to the 
lamina would contain the optic axis of the calcite crystal and the optic axis 
of the twinning layer as well. The plane under reference would thus be a 
plane of symmetry for the entire crystal including the twinning layer. As 
in the case of potassium chlorate recently discussed in these Proceedings,’ 
it can be shown from elementary considerations that when the plane of 
incidence of the light coincides with this symmetry plane, reflections and 
refractions would disappear and the incident light-waves would emerge from 
the twinning layer with their direction of propagation as well as their intensity 
unaltered. 


The foregoing conclusion is readily tested by viewing a distant source 
of light through the calcite rhomb held in front of the eye of the observer 
in such manner that the twinning layer is vertical while the plane of observa- 
tion is horizontal; the source should be viewed through the faces of the 
rhomb which are bounded by edges of which none is parallel to the twin- 
ning layer. Only the undeviated image of the source can then be seen. But 
if the crystal is slightly tilted about a horizontal axis so that the plane of 
incidence of the light ceases to be horizontal, the deviated images make their 
appearance, and rapidly gain in intensity as the plane of incidence is turned 
away from the horizontal in either sense. 

An interesting modification of this experiment is to view an extended 
area of illumination through the faces of the rhomb held in the manner 
indicated above and to rotate the rhomb about a vertical axis so that the 
plane of the twinning layer remains vertical but alters its orientation with 
respect to the direction of observation. A division of the field of observa- 
tion into two parts differing in their intensity of illumination then comes 
into view. A dark and sharply defined boundary is seen to separate the 
two parts of the field except where it is cut by the horizontal symmetry plane. 
This dividing boundary becomes more and more conspicuously observable 
as we move away from the symmetry plane in either direction. Interposing 
a polaroid in front of the rhomb, it is found that the division of the field into 
two parts and the dark boundary separating them are most conspicuous 
when the vibration direction of the polaroid is vertical, but disappear when 
the same is horizontal. It is thereby made evident that the boundary is the 
one indicated by the geometric theory as the limit of the angle of incidence 
of the ordinary waves on the twinning layer beyond which the extraordinary 
reflections and refractions disappear. That the boundary vanishes in the 
symmetry plane and becomes more conspicuous as we move away from it 
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is to be expected, in the circumstances. The boundary is sharply defined 
evidently as a consequence of the weakening of the ordinary transmission 
by reason of the extraordinary reflections and refractions reaching their 
_maximum intensity just prior to their disappearance at the limiting 
incidence. 


A slight modification of the arrangements described above enables us 
to observe the boundaries at which the reflections and refractions of the extra- 
ordinary waves as ordinary waves terminate. For this purpose, the illu- 
minated screen is held to one side so that the light reflected or refracted by 
the twinning layer comes into the field of view. The boundary is seen then 
as an arc of light bordered by coloured fringes and separating a bright from 
a dark part of the field of view. The bright arc is however interrupted at 
the symmetry plane, and gains in intensity as we move away from that plane 
in either direction. By inserting a polaroid before and after the: calcite 
rhomb respectively, it is readily demonstrated that the luminous arc has its 
origin in the trains of extraordinary waves which are transformed to ordinary 
ones by reflection or by refraction. The sharpness and brilliancy of the 
boundary are to be expected alike on geometrical and physical considera- 
tions. 

4. INTERFERENCE PHENOMENA 


As has been already remarked, the deviated images of a distant light 
source seen through the twinning layer are spread out into spectra by reason 
of the difference in the dispersive powers of calcite for the ordinary and extra- 
ordinary waves. Using a narrow slit as the source of light, these spectra 
can be readily observed and are found to be chanelled by interference bands 
crossing them. The position and number of the interferences vary with 
the specimen under observation and thus evidently depend on the thickness 
of the twinning layer in each particular case. Careful examination reveals 
that we have two sets of interferences simultaneously present in the spectra, 
one of which is on a finer scale than the other. The explanation of these 
effects is not far to seek. Returning to the Text-Figs. 1 and 2 above, it 
will be evident that in the passage of the incident waves through the twinning 
layer, the birefringence of the latter would come into play. Hence, each of 
the wave-trains reflected or refracted by the layer would itself be the resultant 
of the wave-trains of the two different species possible within the layer, the 
optical paths of which would be different. Interferences of the same nature 
as those exhibited by crystalline plates in polarised light would therefore 
arise and would manifest themselves in the spectral character of the reflected 
and refracted wave-trains. Effects of the same nature as the ordinary 
interferences of thin plates would also simultaneously appear. The colours 
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observed with extended light sources are a consequence of both species of 
interference. 


Interference phenomena are also strikingly manifested at and near the 
boundaries of transmission, reflection and refraction considered in the 
preceding section. They take the form of alternate dark and bright bands 
of varying sharpness and intensity fringing the boundaries referred to. There 
again, we observe two species of interference, one of which is on a finer 
scale than the other. These boundary effects are best studied using a diffus- 
ing screen illuminated by a monochromatic source, though indeed they can 


be also observed with white light. 
5. . EFFECTS OBSERVED WITH MULTIPLY TWINNED CALCITE 


We now turn to a consideration of the phenomena, typical examples 
of which are reproduced as Figs. 1, 2, 3 and 4 in Plates I and I accompanying 
the paper. The photographs reproduced were obtained using a_ small 
aperture backed by a tungsten filament lamp as the source of light. A photo- 
graphic camera placed about three metres away was fecussed on the source, 
the calcite rhomb being placed immediately in front of the lens. Panchro- 
matic plates were employed for recording the patterns. The brilliant colour 
effects actually observed are, of course, not reproduced. 


Figure | shows the three images cbserved with a calcite rhomb having 
a single twinning layer. {n addition to the bright images, we observe a sys- 
tem of bright streaks, as also some diffuse light surrounding the images. 
The surfaces of the cleavage rhombs had not been specially polished and 
their unevenness was no doubt in part responsible for some of the unwanted 
effects. Nevertheless, our studies with numerous specimens left little doubt 
that the long bright streaks appearing in the photographs and also seen 
visually had their origin within the calcite rhombs. It appears reasonable 
to attribute them to a lack of perfect smoothness of the twinning layers 
within the crystal. 

Another remark may usefully be made with regard to Fig. 1 in the 
Plate. This was recorded through the faces of the calcite rhomb which had 
two of their edges parallel to the twinning layer. In these circumstances, 
the plane of incidence of the light on the twinning layer was, in general, 
far removed from the symmetry plane. By turning the crystal about one 
edge, the images could be made to approach towards or recede from each 
other. Simultaneously, there were large changes in the colour and intensity 
of the images. At particular settings which were nearly but not quite coinci- 
dent, the two outer images vanished but on further rotation they reappeared 
in positions closer still to the principal image. The observations indicated 
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that for the settings at which the images disappeared, the plane of incidence 
was nearly coincident with the symmetry plane. 


Figure 2 in Plate I shows the pattern of a type frequently observed by 
us and which arises from the simultaneous action of twinning layers parallel 
to two of the edges of the rhombohedral cleavage. A noteworthy feature 
in the photograph is that the middle spot in the outer rows has split into two 
components. These were found to be polarised in perpendicular directions. 
Figs. 3 and 4 represent more complicated cases. It is to be noted that the 
spots appear in groups of threes, with the middle one in each row distinctly 
separated into two components. All the patterns exhibited special features 
of polarisation analagous to those observed in the simple case of a single 
twinning layer. 

6. SUMMARY 


Arrays of multi-coloured images of a distant light-source arranged in 
geometric patterns are exhibited by calcite rhombs traversed by twinning 
layers in several directions simultaneously. The paper reproduces four 
photographs of such patterns and also discusses the theory of the reflection 
and refraction of light by twinning layers in calcite. The large difference 
between the ordinary and extra-ordinary dispersive powers results in a linear 
source observed through such a layer appearing drawn out into spectra 
chanelled by interferences. Interference bands also fringe the sharply defined 
boundaries seen separating the different parts of the optical field of reflection 
and refraction. 


7. REFERENCES 


1. Christiaan Huygens .. Treatise on Light, English Translation by Silvanus 
Thompson, University of Chicago Press, 98. 

2. Mascart .. Traite’ D’Optique, 2, 

3. Grailich .. Pogg. Ann., 1856, 98, 205. 

4. Osthoff .. N. Jahrb. f. Miner. Beil-Bd. 20, e 

5. Pockels .. Lehrbuch Der Kristalloptik, 203. 

6. Rayleigh .. Scientific Papers, 3, 212 


7. Raman and Krishnamurti... Proc. Ind. Acad. Sci., 1953, 38 A, 263. 








ON THE DISTRIBUTION OF THERMAL 
NEUTRONS FROM A FAST NEUTRON SOURCE 
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By R. RAMANNA* AND S. B. D. IYENGAR 
(Tata Institute of Fundamental Research, Bombay) 
Received April 23, 1954 
(Communicated by Dr. H. J. Bhabha, F.R.s., F.A.SC.) 


1. INTRODUCTORY 


THE object of the investigation is to determine the effect of a paraffin reflector 
on the spatial distribution of thermal neutrons due to a point source of fast 
neutrons kept at the surface of a graphite moderator. If the effect of the 
reflector is to reduce the leakage of neutrons from the sides of the graphite 
without appreciably reducing the total flux due to its own absorption, it 


may be possible to avoid the use of expensive graphite for reflector and use 
paraffin or water instead. 


In order to determine the effect of the reflector, the thermal neutron 
density in a graphite medium due to a point fast neutron source placed on 
its surface and covered by (1) a paraffin reflector and (2) a graphite reflector 
was measured. An expression was obtained for the thermal neutron flux 
in a medium due to a point fast neutron source on the interface between 
it and another medium, by means of the two group theory. The experi- 
mental results were compared with the theoretical distribution. The results 
were also compared with a theoretical distribution for a uniform plane fast 
neutron source, by integrating the results for the point source at various 
points. The agreement between theory and experiment is discussed. 


2. THEORY 


Let A be a cylindrical block of moderator of radius R and length L, 
and B a cylindrical block of reflector of same radius and length L,. At 
the centre of the interface is kept a point source of fast neutrons of strength 
S,. If we suppose that the neutrons in the medium A are divided into 2 
groups, one fast and one thermal, then we can write for the fast group flux, 


Diy bra hak aA $:,=9, (1) 


* Reactor group, Atomic Energy Commission, 
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where the subscript 1 refers to the fast group, 
¢,, 1s the fast group flux, 
=, the slowing down cross-section, and 
D,, the fast diffusion coefficient. 

For the thermal group flux, 
Do, Vv bo, — Xn Faq + 21, $1, = 9 

where the subscript 2 refers to the thermal group, 

do, is the thermal group flux, 


2, is the macroscopic absorption cross-section 
and 
D,, the thermal diffusion coefficient. 


Similarly, in the medium B, 
Diy Voip wa ~'1p Pin = 0, (3) 
Dox v" Pop = aon Poy =e Jip Pip, (4) 


where D,,», Dap, 2;, and 2,, are the corresponding constants for the 
medium B. If we write the solution of equation (1) as Z,, we can write a 
solution of equation (2) as Z, +S, Z, where Z, is a solution of 


Don V2, — Zea Goa =0. 

Similarly we can write the solutions of (3) and (4) as 
Z, and Z, + S, Z,. 

By substitution we have 


Py 1 
S cmead 2% 9 9 3 
" Dox Per? — Pra” (5) 
> ] 


Ss, = wat | 3 i. 
° Dey, Poy” — Pir” (6) 


= x14 = Joa _ a1 = aon 
Pia D.’ Pua > |e Pia. > Nhs Pip >= Do» (7) 


where 


At first we shall suppose that the source is not a point source but a plane 
source distributed over the interface, the strength of the source being a func- 
tion only of the distance from the centre. Let S (r) be the strength per unit 





10 R. RAMANNA AND S. B. D. IYENGAR 


area of the source at distance r from the centre. The equations are all 
essentially of the type (1), and are solved by separation of variables in cylin- 














MEDIUM A 
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drical co-ordinates taken as shown in Fig. 1. Due to symmetry about the 
z-axis, there are only two variables. We obtain 


$1, = (Ay, e*? + A’, e-%] [Cyn Jy (kr) + Cin Yo (kr)] 
do, = [Age + A’ay e777] [Caz J, (kr) + C’o. Y, (kr)] + S, $1, 


and similar expressions for ¢,,, and 4, where 


a= Vk? + p,.2, B= Vk? + pips y= VA KE + Po”, 8= Vk? + pox? (8) 
the solutions being summed over all possible values of the separation 
constant k. Our solutions must satisfy the following boundary con- 
ditions: 
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(1) The ¢’s should be finite at all points. The coefficients of the Y, (kr) 
terms should therefore be zero. 


CRr= Cn = 6, ete. 
(2) At z=Ly,, ¢,,=0, ¢., =90, for all values of r, i.e., 
Aye — A’ 072", Asp = — A’op e72™ 
At z=— Lg, ¢:,=9, ¢.,=0 for all values of r, i.e., 
Bs. = — By, e727, B’s;, = — Boy e778 
(3) At r=R, ¢,,=¢2,= ¢1—3= de3= 0, for all values of z, i.e., 
J, (AR) =0. 


k can only take values k,, ky, .... 


where these are the solutions of above equation for k. Making use of the 
above relations our solutions can be written as 


by, = S Ayn [eman® — e2n(2-2L4)] J, (Kr), (9 a) 


n=1 


co 


Por =F Asn - or ieee, J, (kyr) 2 2 S, py As (9 b) 


n=1 


co 
y’ 


dip= XY Byn [eFn® — e-Palz+ 2nd] J, (kyr), (9 c) 


n=1 


gon= F Ban [e5"* — e-ul2+210)] J, (kar) + Sp di. (9 d) 


=1 


(4) At the interface, i.e., at z=0 


3 


$4, =p, and do,—= doy for all values of r. Hence 


E Ain (1 — e-®anls) J, (Kknr) = Bin (1 — e-2R aL) Jy Knr), 


> Asn (1 — e-*ha) J, (Kn) +S, > Aw a= e-*anla) J, (nr) 


n n=1 


n=1 


= F Bon (1 — e776") Jy (Knr) + Sp 5 By (1 — e-*Bnln) Jy Knr). 


Multiplying the first relation on both sides by r J, (kmr)eand integrating 
with respect to + between 0 and R, we get 
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co R 
D Ain (1 — e-*@eh) fr I Kmr)So (Knr) ar 


n=1 


a > Bin (1 — e7? Pal) fr Jy (Kmr) Jo (kyr) dr, 


which simplifies to 
Aum (1 — 724ml) = Bim (1 — e7? Pmt») 
Similarly from the other relation 
Aom (1 — e7?%mE4) + S, Aim (1 — e-2amls) 
= Bom (1 — e7-5™™) + §, Bum (1 — e72P mix) (11) 
(5) At the interface, i.e., at z=0, 


Ghz — Giyz=S(@), and (Io.)2 — Jey)z=9 for all values of r, (J,,), being 
the z-component of the current density J,,, etc. 


d 
Siadz ai Dia Ps » etc. 


From (9), we then have 
Siaz= Dig 5 On Ary [e-*"* + e%'2-214)] J, (kyr). 


We obtain similar expressions for (Jo,)2, (Jzp)2 and (Js»)-. 
We then otain 


Di, z Qn Ayn (1 + e-26nb) J, (Kur) + Dip > Bn Bry (1 + e-7Fal) J, (Kar) 


n=1 n=1 


= §(r). 
Deon y Yn Aon (1 +e? 74) Jy (Knr)+Dea Sa E Ain G + a ™")5, Gy) 


n=1 n=1 


Dos > bn Bon (1 +-e~*SeL2) J, (kyr)—Dep Sp E Bn Bin (1 + e7Fa!4) J, (Knr) 


n=1 


Multiplying both sides by r J, (Kmr) and integrating with respect to r between 
0 and R and simplifying, we obtain 


Dis am Aim (1 cae e~2amL«) oe Dip Bm Bim (1 -- e~2BmLr) 


1 Sr) rJy Kmr) ar 
~ b km? R? (km RP 
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Dea ¥m Am (1 + e?%™) + Dox Sam Arm (1 + e-*am™s) 
= — Day 8m Bem (1 + e775") — Dag Sp Bm Bim (1+e-7F mt) 
(13) 


Solving (10), (11), (12) and (13) for the constants and substituting in (9) 
we get 


oo 


bis = DT aaa [eH — 8 e0)] Tyr) (14a) 


n=1 


co 
Pin = eas [ef — ert 2+ 01) Jy Cnr) (14 b) 


n= 


r= 1 
co 

oo, = c. [ An (e- In — ey (2—2La)) 
n=1 


S, (e-o a et (Z—2La)) 


1 — e-2an La | Jo (Knv) 


co 
po, = > Cc, [Bn (e5n2 — e-dn (2+2L0)) 
n=1 


S,, (ent — e-Be (2+2Ls) 


+ SEN papas — | Jo nr) (14 d) 


tie i 2 FS (r) Iq (kyr) dr 
nn SS eee RAO : 
[D,, 4» coth (ay Ly) + D,p8y coth (Brly)n?R7[5, (k,R)}? 
(15) 


i S;Do, [5,, coth (S,Ly)—Bp coth (BnLy)J—S, [D2.¢n coth anLy+ D,.35n, Coth (8,L,)] 
’ (1 — e-?%'4) [Dey yn COth (Yn La) + Day dn coth (5, Ly) 











(16) 


With a similar expression for B.,,. 


To apply these results to the case of a point source, S (r) must be of such a 
form that 


S (r)=0 when r + 0 


and the total number of neutrons emitted by the interface should be equal 
to S, 


[S(r):2ardr=S, 
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if 3$@= Se! = , O(r) being the Dirac Delta-function, these conditions are 


satisfied. Substituting this expression for S (r) in equation (15), we obtain 
a 
7 [D,, 4n Coth (a, L,) + Diy By coth (By, Ly)] ky? R* [Jy (AnR)/? 
With this value of C,, equations (14) give the final solution. 
The problem was also solved for a uniform plane source distribution 


on the interface. There is now only one co-ordinate z and the cylinders 
become of infinite radius. Let S be the strength of source per unit area. 


The solutions are: 


‘ 


1, = =a [e-Psz oo eP1a(2—2La)] (18) 


S, (e-Pz — eP1A(2—2La)) 4 


de,=C | A (e~ P= — eP2a (2-2La)) + 1 — e-2PiaLa jo 


where 
Dy, Pra COth (Py, Ly) + Dyy Pyy Coth (py, Ly) 
ei mn (S,—S,) Dey P2,coth(p2 Lb J—[DorPra coth(p,,L4)+Deppiz, coth (PiyLy)] 
(1 — e-*Phs) [Do, Poy COth (Po,L,) + Dep Poy COth (Poy Ly)] 
EXPERIMENT 
The experimental arrangement for the measurement of the thermal 
neutron distributions is shown in Fig. 2. It consisted of a large graphite 
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block 120 cm. wide, 80cm. long and 47 cm. deep with a rectangular cavity 
of depth 30cm. from one surface for introducing the detectors. The space 
was then refilled with graphite bricks and covered by a paraffin or graphite 
reflector as the case may be. The neutron densities were measured by means 
of the activation of Mn-Ni foils previously calibrated for relative efficiency. 
The foils were placed one below the other inside slits cut in a graphite brick 
and then placed in the required position inside the rectangular space. A 
Ra-Be neutron source of 50 mec. strength was placed in various positions 
along the reflector-graphite interface and the activation of the foils for each 
position of the source was measured. The counting assembly consisted of 
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a pair of geiger counters in parallel shielded for cosmic rays and a conven- 
tional counting equipment. 


It will be seen that the geometry of the experiment does not exactly 
coincide with the geometry assumed for deriving the theoretical expression, 
The graphite and paraffin blocks are not cylindrical, the source is not always 
at the geometrical centre of the interface and the detectors are also not along 
the central axis. But due to the large size of blocks, the exact geometry 
ceases to be of importance. The source and detectors though asymmetrically 
placed are not far from the central axis, there being at least 35 cm. of graphite 
surrounding them on the thinnest side. So the theoretical expression may 
be applied to explain the results of the experiment. 
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Fig. 3 shows the variation of Log. (thermal neutron flux) with z, 
for various values of r. 


The experiment was also performed with a graphite reflector but for only 
one position of source with r= 0. The Log. (thermal neutron flux) is plotted 
against z, for r=0, for both graphite and paraffin reflector in Fig. 4. The 
paraffin value is higher than the graphite value for small values of z, but at 
about z= 15-5 cm. it becomes equal and then becomes less than the graphite 
value. 


In order to determine the flux distribution for a plane source, the sur- 
face around the point r=0, was divided into concentric rings, of width 
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10cm. corresponding to the various values of r, and the contributions to 
any point along the z-axis from these various rings were integrated and the 
flux due to a plane source disc of 55cm. radius was obtained. Fig, 5 
shows the Log. (thermal neutron density) points so obtained and also the 
theoretical curves obtained from expression (19) for paraffin. The agree- 
ment, it will be seen, is quite good. 


Values of ¢,, the thermal neutron density for various values of r and 
z were calculated from expression (14c). The constants for paraffin and 
graphite used for the purpose are shown in the appendix. 
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Fig. 6 shows a polar diagram of Log. ¢,, plotted against 0, the angle 
between the line joining the point to the source and the z-axis, for points 
lying along straight lines parallel to the z-axis at various distances r from it. 
The diagram shows theoretical curves for r—0, 10, 20, 30, 40 and 50cm., 
with experimental points also indicated. It can be noted that for r=0 the 
curve is just a straight line. It will be seen that the agreement between theory 
and experiment is good for large values of z and r but poor for small values. 
This may be explained as due to the fact that the assumptions on which 
the theoretical result is derived, such as that simple diffusion theory applies 
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throughout the region and that the source is a point source, become unten- 
able at small distances from source. 


A calculation from (14 c) of the flux for the paraffin and the graphite 
reflector for ry =0 showed that the graphite values are always much higher 
than the paraffin values. This led us to consider that the graphite used may 
not have been quite pure. The values of the constants mentioned in the 
Appendix, are for graphite with impurity less than 2 parts per million, the 
impurity being wholly due to boron content. If greater contamination of 
boron is considered, the other constants remain practically unchanged but 
the value of p, increases. Calculations were made for various values of 
p, and Fig. 7 shows the theoretical curves for a paraffin reflector for three 
values of r and for a graphite reflector for r=0, for an arbitrary value of 
Pp, =1. It isseen that these curves agree qualitatively with those of the 
experimental curves given in Fig. 4 only for this very high value of pg, the 
paraffin value being more than that of graphite for lower value of z and 
less than the graphite for higher value of z. A closer fit to the experiment 
was not attempted as the exact value of p, to fit the experimental results 
was not of any particular importance. 


It appears that on the two group theory for pure graphite the use of 
the paraffin instead of the graphite reflector reduces substantially the thermal 
neutron flux, the intensity at z= 0 for the paraffin reflector being only 0-34 


times the value for the graphite reflector. 


For larger values of z this ratio 
continues to have almost the same value. 


The experimental part of this work was done at the Centre d’Etudes 
Nucléaires, Saclay, Paris. One of us (R.R.) has pleasure in thanking the 
Commissariat a |’Energie Atomique, France, for providing the facilities to 
perform the experiment and also in thanking M. Martelly for many useful 
discussions. Our thanks are also due to Dr. H. J. Bhabha for his 
encouragement during the course of this work. 
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APPENDIX 
Constant 
Slowing down Cross-Section 
Fast Diffusion Coefficient 
Macroscopic abs. Cross-Section . 
Diffusion Coefficient (thermal) 
> 


cat 


A 


Graphite 
0-00302 
1-057 
0-00036 
0-903 


0-0535 


0-01997 ‘) 


—0-37(t) 


+ For impure Graphite these values are different. 


Paraffin 
0-0505 
1-515 
0-017 
0-142 


0-175 


0-346 


0-398 
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HETEROCYCLIC COMPOUNDS* 


Part If. Synthesis of 4-Quinazolones with Substituents 
in the 2- and 3-positions 


By M. R. SUBBARAM 
(University of Madras, Madras—25) 


Received January 2!, 1954 
(Communicated by Dr. K. N. Menon, F.A.sc.) 


METHODS previously employed for the preparation of 4-quinazolones include, 
among others, the treatment of N-acylanthranilides to elevated temperatures! 
or treatment of acylanthranils with primary amines.2, McKee, McKee and 
Bost? employed the action of substituted benzonitriles with acid chlorides 
to get 2-substituted-4-quinazolones. 2: 3-substituted-4-quinazolones were 
obtained by Grimmel, Guenther and Morgan‘ by a method which allows 
the preparation of these compounds in one simple operation at moderate 
temperatures. When primary amines are condensed with N-acetylanthranilic 
acid in toluene or any other suitable solvent, with the aid of phosphorus 
trichloride, 4-quinazolones are formed in fairly good yield. 


As part of a scheme of study of Heterocyclic Compounds in progress 
in this laboratory, the synthesis of some 4-quinazolones (see experimental 
section) was taken up. While this work was nearing completion, Rani, 
Vig, Gupta and Narang’ published a paper on similar lines, embodying the 
synthesis of several 2: 3-substituted-4-quinazolones. The 2: 3-substituted- 
4-quinazolones reported in the present work have a methyl group in posi- 
tion 2, and a substituted aryl group in position 3. The N-acetylanthranilic 
acid required for this work was most satisfactorily obtained by the method 
of Bogert and Gotthelf.® 


EXPERIMENTAL 


2-Methyl-3-(2'-chloro-4'-nitrophenyl)-4-quinazolone—In a 500 ml. three- 
necked flask equipped with stirrer, reflux condenser and dropping funnel, 
and mounted in an oil-bath, were placed 1-8 g. (0-01 mole) of N-acetyl- 
anthranilic acid, 1-7 g. (0-01 mole) of 2-chloro-4-nitroaniline and 75 ml. 
* Part I. Menon and Raman, Proczedings, 1953, 38A, 132. 
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of dry toluene. The mixture was stirred and treated drop by drop with a 
solution of 0-5 g. (0-0034 mole) of phosphorus trichloride in 5 ml. of toluene 
over a period of five minutes. The resulting suspension was stirred and 
refluxed for two hours. After cooling, the contents of the flask were trans- 
ferred to a round-bottomed flask, treated with 100ml. of a 10% sodium 
carbonate solution and steam distilled, to remove toluene. The residue was 
collected and crystallised from acetone giving 1-7g. (63% yield) of the 


quinazolone, m.p. 170-71°. (C,;H,,O,;N;Cl: Calculated—N, 13-31. Found: 
13-35). 


2-Methyl-3-(2'-methyl-5'-chloropheny!)-4-quinazolone.—Condensation of 
1-8 g. (0-01 mole) of N-acetylanthranilic acid with 1-4g. (0-01 mole) of 
2-methyl-5-chloroaniline gave 2-0 g. of the pure quinazolone (from acetone), 
m.p. 142-44°. (C,,H,,ON.CI: Calculated—N, 9-84. Found: 9-59). 


2-Methyl-3-(2’ :6’-dimethylphenyl)-4-quinazolone.—N-acetylanthranilic acid 
(1-8 g., 0-01 mole) and 2: 6-dimethyl aniline (1-2 g., 0-Olmole) gave 1-8 g. 
of the pure quinazolone (from acetone), m.p. 130-32° (C,;H,,ON,: Calcu- 
lated—N, 10-60. Found: 10-54). 


2: 6-Dimethyl-3-(0-toly!)-4-quinazolone.—Molar proportions of N-acetyl- 
5-methyl anthranilic acid’ (2-0 g.) and o-toluidine (1-0 g.) gave 1-0 g. of the 
pure quinazolone (from acetone), m.p. 159-61° (C,;H,,ON,: Calculated— 
N, 10°60. Found: 10-59). 


2: 6-Dimethyl-3-(o-anisyl)-4-quinazolone.—Molar proportions of N-acetyl- 
5-methyl anthranilic acid (2-0 g.) and o-anisidine (1-2 g.) gave 1-8 g. of the 
pure quinazolone (from acetone), m.p. 176-78° (Calculated for C,;H,,O.Nz: 
N, 10-00. Found: 10-18). 


2: 6-Dimethyl-3-(0-phenetyl)-4-quinazolone.—Molar proportions of N- 
acetyl-S-methyl anthranilic acid (2:0g.) and o-phenetidine (1-3 g.) gave 
2-1 g. of the pure quinazolone (from acetone), m.p. 138-40°. (C,.H,,ON,: 
Calculated—N, 9-53. Found: 9-53). 


I wish to take this opportunity to express my thanks to Dr. K. N. Menon 
for his interest and guidance during the course of this work. 


SUMMARY 


The synthesis of 2-methyl-3-aryl-4-quinazolones is reported. 





M. R. SUBBARAM 


Korner 
Bogert and Chambers 
——— and May 
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GENERALIZED SINGULAR POINTS WITH 
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SUMMARY 


Generalized singular points are used to discuss irrotational and 
viscous flows produced in an infinite liquid by a moving solid. It is found 
that the irrotational motion of translation is the same as that due to a 
generalized doublet and that of rotation the same as that due to a “* rota- 
tion singular point’. The corresponding viscous problems are solved 
by superposing on the irrotational motion a solution due to a concentrated 
force or a couple. 


1. INTRODUCTION 


SINGULAR points play an important part in all branches of mathematics. 
They appear as sources, sinks, doublets, rectilinear vortices, electric charges, 
magnetic particles, electric currents, etc. It is customary to think of them 
as the limiting cases of spherical or circular elements. It is well known 
that all potential flow can be produced by a suitable distribution of them. 
But, excepting for simple boundaries like planes, spheres or circles, this 
distribution is infinite in number, and, hence for practical applications it is 
not very suitable. It is not generally appreciated that a generalization of 
these points as ultimate forms of a family of closed surface bodies can simplify 
a large number of boundary problems by reducing the infinite number of 
spherical or circular points to a finite number of generalized singular points. 
For example, an ellipsoidal source or doublet can play the same part for 
ellipsoidal boundaries as an ordinary source or doublet does for spherical 
boundaries. In fact the uniform non-viscous flow due to the motion of an 
ellipsoid can be produced by an ellipsoidal doublet. 


It is proposed to discuss how such generalized singular points may be 
obtained and used for flow problems. It is found that :— 


(i) the irrotational motion of a solid through a non-viscous liquid is 
the same as that due to a generalized doublet; 


* Presented at the Eighth International Congress on Theoretical and Applied Mechanics 
held in Istanbul, 1952. 
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(ii) the irrotational motion of rotation is the same as that due to “ag 
rotation singular point ’’; 


(iii) the slow viscous motion of the solid can be obtained by superposing 
on an irrotational motion due to a generalized doublet a solution due to 4 
concentrated force in the direction of motion; 


(iv) the slow viscous motion due to rotation of the solid, if stable, can 
be obtained by superposing on the motion due to “‘ a rotation singular point,” 
the solution due to ‘‘a centre of rotation”. 


The corresponding analogies in electricity and magnetism and heat 
conduction can be easily written down. 





















2. IRROTATIONAL FLOw 
Let us take a family of closed surfaces given by 
f(x, y, 2, &)=0, (1) 


and let us assume that the surface of the moving solid A is given by ¢=,. 


The family given by =a constant, is not necessarily a family of equipotential 
surfaces. In order that it be so we should have 


V? [F(é)]=9, 









which gives that! 


v(x) +G)+G] 


should be a function of € only, say w’ (é). 
the corresponding potential is found to be 


cfexp[—#(é)])dé+d 
c and d being constants. 


When this condition is satisfied 
(3) 


Let V be the potential at an external point due to the solid A whose 
mass is M. Then, by Gauss’ theorem the normal flux across any surface 


enclosing S is 
ff Pai ds = — 4n M, 
on 


the element of normal 6” being drawn outwards. The limiting form of A 
when it reduces to a point or a line will be called the generalized source of 


strength M. Its potential will be taken as V. V satisfies the following 
conditions :— 
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(i) V2 V=0 throughout all space unoccupied by matter; 
(ii) the flow across any closed surface containing it is constant; 


(iii) Grad. V vanishes at infinity. The generalized source thus satisfies 
all the requisite conditions. 


A doublet in the x-direction is obtained by displacing the source or solid 
through#a small distance in the x-direction. The corresponding potential 
isp )V/dx, » being the strength of the doublet. Again, the composite body 
whose density at any point Q is the difference of the densities at Q of the 
given body in two neighbouring positions is equivalent to a boundary layer 
of density p cos @ placed on the surface of the solid A, 6 being the angle which 
the normal makes with the x-axis. Thus the doublet can be taken to be 
equivalent to a boundary layer on the surface of A of surface density 
pcos 8. 


If, instead of giving a translatory displacement, we turn the body round 
the z-axis through a small angle, we get a “ rotation singular point ’’, whose 
potential 


b=n(x dV binds *) . (4) 


oy ox 


is the same as that of a boundry layer on the surface of A of surface density 
u(mx — ly), 1, m, n being the direction cosines of the normal to the surface 
S of the body A. 


We shall now show how these singular points may be used to discuss 
the motion due to the translation and rotation of the solid A in an infinite 
non-viscous liquid. 


At first we show that a generalized doublet or boundary layer gives 
the flow due to the motion of the solid A through an infinite non-viscous 
liquid. 

By Green’s theorem, if N, N’ are the normal intensities due to the boun- 
dary layer of density p cos 4 at points just inside and just outside the layer, 
we have 

N’ — N= 4rp cos 6. (5) 


Now, since the surface S is an equipotential, we can take N=0. Thus we 
see that the normal component of the potential due to the doublet assumes 
a value proportional to cos @, which is the condition to be satisfied by the 
potential in the corresponding irrotational motion. The constants can be 
easily adjusted in terms of the velocity of the solid A. 
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We have thus shown that the motion of the solid A through an infinite 
non-viscous liquid can be represented by a generalised doublet, whose poten- 
tial can be obtained from the static potential due to the solid A. 


In the case of a sphere moving through an infinite liquid we have the 
following results :— 


The potential V at any external point is 


M 


r 


The potential due to the corresponding doublet is 


a d 5 oo: 
ees 


which is the known result for the potential of liquid flow. 
In like manner for a circular cylinder we have 


V=c— 2Mlogr, 
so that 


px 
o="Y, 
which is again a known result. 


Let us take the case of an ellipsoid. The confocal family of ellipsoids 
is given by 


The potential due to an ellipsoidal source is? 


co 
- ; 
v= (1 A ee du 
atu B+u cosy) 4° 
K 


@ 
The potential of the ellipsoidal doublet is therefore given by 


dV 
¢ =—# = =HXe, 


#,, being a constant, and 4, a are given by 


= (a? — r du 
A=(a? +0 +n (2+), f arya 
& 


which result is again known. 
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For an elliptic cylinder we have* 


oe. fF | 
re ae - Da <6, a (+. Be a : 
V oP ax rae 7) mpab log (a’ + b’) + « (7.1) 
a’, b’ being the axes of the confocal ellipse through the attracted point. 
This value of V gives 


o=p —, ab of , a =c cosh é, b’=c sinh &, 


a+b’ a 
x=ccosh écos 7, y=c sinh € sin 7. 
Substituting these values we get the known result 
=p, exp (— £) cos 7. 
For the elliptic source V is of the form* 
V = — 4 xpab (e—*é cos 2n + 


(7.3) 
For the motion of rotation, we easily get the following known results 
for the value of ¢ given in (4) :— 
(i) for a sphere or a circular cylinder ¢=0; 
(ii) for an ellipsoid* ¢=p, (8 — a) xy, 
(iii) for an elliptic cylinder ¢=, e~*é sin 27. (8) 
Thus we see that the irrotational motion of an infinite liquid due to the 


translation or rotation of the solid A can be easily deduced from the static 


potential of the body at any external point, by using the idea of generalised 
singular points. 


For any internal point V satisfies Poisson’s equation 
V2 V=-— 4a or — 2p (9) 


according as it is a three-dimensional or two-dimensional problem. p is assumed 
to be a constant, and hence V does not satisfy Laplace’s equation. But its 
first partial derivatives with respect to x, y, z satisfy it. The potential for 
the transational and rotational flows of the liquid comment in the surface 
S of the body A remain of the form 


dV dV dV 
t= and p (x ~~ y =) , (10) 


* g is obtained from a in (6.3) by interchanging @ and b. 


A ee 
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The following results now follow immediately :— denot 
(1) for the motion of translation along the x-axis of any solid 6 = — Ux surfa 
(2) for the rotation of sphere or circular cylinder 6=0; ! 
(3) for rotation of an ellipsoid, we have to use V=D— Ax? — By? we g 
— Cz?, D, A, B, C being constants. This gives 
d= pyxy 
(4) for an elliptic cylinder Sine 
= ab (x* , y? Pe 
7 “Patb\a be ob = mXy. we | 
Thus liquid motion contained in a closed surface given by (1) can be 
deduced from the static potential at any point inside the body A. 
3. Viscous FLow 
u,v 
For the slow motion of a solid through a viscous liquid we have to 
consider generalised singular points of various types. In the first type we 
determine the displacements produced in an infinite elastic solid wherein 
a concentrated force X, acts at the origin in the direction of the X-axis. In 
this case we know that the displacements components can be put in the 
form® The 


“= 


RY) dH IG 
dx + dy —_ 2’ etc. (11.1) 


and the body force X, Y, Z in the form 










x= 


where 





(A+ 2n) V°d + p® =0, ny*F + pL=0, »y*G + pM=0, 


- et 
ey*H + pN=0. (12) 


th 
The Laplace’s equation in orthogonal curvilinear co-ordinates £, y, é is 


d 1 OV — an 
Y; (hats = 2) + two similar terms = 0, 


1 (x, Pry, fr2V 
hy? = (3) + (3 | (32 — (13) 
Let the surface of the moving solid A be given by the equipotential 
= &,. There is a solution of (13) which is a function of € only. Let it be 


where 
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denoted by %#(£). It should be such, that its flux across any equi-potential 
surface is equal to 47 in three dimensions and 27 in two dimensions. 


Putting 
p J SJ X’ dx’dy'dz' = X,, 

we get? 
| 


G=—_ X 
p 


: ~ int wa! x * Be 
7T 


0 dx 4p 0 IZ 


Since a 
Vv" (xt) = 2 a 
we can put 


8 (2a) xe, P= =@, G=-— Xo 


= 


d 
_— Sar 


_ X , 
H= Sr ye. 


u, v, w now take the form 


(A a L) Xo Pi) 
Sm (A + 2) dx 


(A+p)Xo (2 r) 
Sap (A+ 2") \dp’ dz 


X 
“= — v0 


(xp) + 


Amp 


) (xz). 


w= 


These give 
uu, WwW, WwW Xo de 
x «(Oy Oz ~ Oe (A + + 2p) ax ? 
— (3A + 4u) Xy rp _ (A + #) Xo >? (xa) 
7 47 (A+ 24) rx 4m (A + 2p) dx? * 
— (A + #) Xo a? 
"ye = ~ a (X+ Qu) dyaz OH) ete. — (17.2) 


The equations of equilibrium are satisfied with these values of ryx, Ty, 
etc. They have to satisfy the condition that across any surface é = constant, 
they give a force X, in the x direction, independent of é. 


X,, the traction component along the x-axis is 
Xy = [tyx + Mtxy + NT xe 
a (A + p) Xy o2 >? 
4a (X + 2p) |! ax? T ™ yyy 12 sz () 


Xo d | o* | A+ p dep 
v dn (755 ) ns, )¥+ Qa (A + 2p) 0 ax ° 


(18.1) 
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/, m, n being the direction cosines of the normal drawn to &= constant. 
If Sy is an element of this normal we can write (18.1) as 


_ (A+ bt) Xo . do / oe se uXy deb 
? 4a (A + 2)” dy ) © 4 (A +2p) dy 


(A + LL) Xo / yp 
4n(A+2n) dx° 


We shall now show that 
ff X,ds = Xo; 
the integration being taken over any closed surface €= constant. 


We shall assume that the liquid is bounded by the surface S (é=é,) 
and an equipotential surface at infinity given by €= oo. In the usual notation 


ff X,ds = ff X,ds— ff X,ds. (19.1) 


But X, is given by (18.2). Also 


Jf- » Gre Sf x; a» (Ge) 


- ffi . ds — {fi =: d 


But, since 4-»+1/r, at infinity, we have 


n) (*) a 7 
ise a | x) a = 37 {fi ee oe 
ff-2) (=) aes ffrx ie & be. 


ffx oF ds — (f wr ds =0, 
Sf ae ds = — 4n. 


Similarly 


so that 
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Thus 


fy yds = — on a Sa = — X,, (19.5) 


which shows that a force X, should be applied in the x-direction. Similarly 


y,=— ATH Xo x ( (A + #) Xo , 
y 4a(A + 2p)” dy v has —" —< 


(A + KL) Xo _ 9 (*¥ ie (A + bh) Xo 1 def (19.6) 


me" (A+ 2u) ~ dy\ az) 4r(A + 2p)” ox 


from which we easily see that 


SS Y,ds=0, ffZ,ds=0. 


For the vanishing of the couple we should have 


SJ(zY, — yz,) ds=0, (20) 


ose (/z + nx) — = (ly + mx) — ny “ oe 2 = ds =0. 


ox 


Transforming into volume integrals, we see that this condition is also 
satisfied. 


For a liquid A~+0, A-+oo, such that A A—-— p. Thus we get from (16) 
0? 0 
—_s fe (x) — 2%), 


Tyy= — P+ 2B (x), 


> 


>" 
T2z= —p+ 2B yz? (xi), 


2 


re >? 
rye= 2H B soe (x9), Tor = 2B >, (ay) — WB ~ 1% 


~ (xy) — 2B OY (21.2) 


Txy= 2 B oy 


ap 


Where y*s=0 and B= — X,/8mp. 
A3 
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The tractions becom 


The displacement u, v, w take the form 
d) 
_ : “ ) 
u=B [ = (xes) — 2b |. 


t= (: 


oN 
dy’ WZ 


) (xe). (21.4) 


Let V be the potential of the corresponding irrotational flow. On the 
boundary of the moving solid (= constant) we have 
dx dV dx 


dE? or U 22) 


dV 
—h, dé =h,U x ~~ U 3: 22 


if V is of the form xa (€), the velocities for the slow viscous motion can be 
at once obtained by superposing on the generalized doublet of the correspond- 
ing irrotational motion a generalized singular point of the first type. Thus 
we get 


AN XV =X d «4 
me dx Bam E (2) 24 ; 


WX 9 
oy Sm dy 


XV =X d 
dZ Simp dZ 


vy =A (x), 


w=A (x), 


where the constant A can be determined from the boundary conditions. 


If the irrotational motion is represented by harmonics of the type Vj, 
V in (23) will have V, and all harmonies of an order less than V,. 


When V= xa (&) we get from (22) 


A [ay + xa) 6] — gh0 | ye’ 3 — yy | =U 
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where ao, %, are the values of a, 5 on the surface of the solid, dashes denote 
differentiation with respect to é, and U is the uniform velocity of the solid A. 
Thus 


Aag+ we =U, Aa’= aa ; (25.1) 


These values of X, and A also ensure that v=0, w=0O on the boundary. 
The drag suffered by the solid and the constant A are given by 


— 8p ay U i Sn bo , (25.2) 
" Pag + fg a, Boag + pq ao Sie 


The second case when V is to be represented by a series of harmonics can 
be similarly treated. 
In the particular case of the sphere we get 


F 2x I a ' 
V=x = 3 pa? b= és _ a = g 2°, X= bray U. 
These give the known values of u, v, w 


For a circular cylinder we get® 


x 
View we: me, * = — logr, 


Amp a9’U 8m U 
Pe ee a er ; (27) 
Poty + ty ao l 2 loga 


An infinite constant is neglected in the value of %, which is a well-known 
drawback. 


For an elliptic cylinder we have 
a=e-€ sech é, 


x, 4m U(a + 6) 
° a—€&j(a+ bd)’ 


a=c cosh &, b=c sinh €,, c?=a? — b?, 
For an ellipsoid we have the following results :— 


. l : l Ua? 
fy =— »% =—-y , A= = ae 
2abc a®bhe 29 + aa, ’ 


— 16mU 
X= 2%, + a*a, : 


and the values for u, v, w which are known.’ 
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For the slow rotatory motion of a sphere or a cylinder in a viscous liquid 
we have to use the generalised singular point, called “* the centre of rota- 
tion’. In both cases there is no irrotational motion, and hence we can take 
V=0. For a sphere, if G is the couple of rotation, the corresponding dis- 
placements in an incompressible infinite solid are 


u, v= es ( . — =) : ,w=0 (30.1) 
If w, is the angular velocity and a the radius of the sphere 
G = 8m a’ a, (30. 
For a circular cylinder of radius a we get the following results :— 


G 


u,v = 
, Amp 


3 3 oe 
of =e) logr, w =0 (31. 
G = 4a? a, (St. 
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ON THE SEMI-EMPIRICAL MASS-FORMULA AND 
ALPHA-DISINTEGRATION ENERGY IN THE 
REGION OF HEAVY AND MEDIUM 
HEAVY NUCLIDES 


By G. P. Dupe, F.A.Sc. AND LAL SAHEB SINGH 


(Department of Physics, Patna University, Patna) 


ABSTRACT 


In the present paper it has been shown that Fermi’s semi-empirical 
atomic mass-formula is inadequate for the estimation of alpha-disintegra- 
tion energy in the medium heavy and heavy nuclide-regions. A suitable 
correction-term has been added to the mass-formula based on Duckworth’s 
new, atomic masses and similar to the one proposed by Stern for heavy 
elements (A > 208). It is seen that the calculated alpha-decay energies 
are in numerical agreement with the observed data except for the magic 
number nuclei which can be satisfactorily explained on the nuclear shell- 
structure. But even in the latter case the agreement with the experi- 
mental results is better when the correction-term is used. 


INTRODUCTION 


Tue determination of atomic masses to a high degree of precision is of great 
importance because of its many significant applications in Nuclear Physics. 
The binding energy of the nuclei, conditions of nuclear stability, radioactive 
properties of isotopes, excited states of parent and daughter-nuclei, their 
inner structural informations, the nature of nuclear forces—all these need 
directly or indirectly for their estimation sufficiently accurate atomic masses. 
The discrepancy discovered in any one of them requires significant inter- 
pretations. The accuracy obtainable in the most refined mass spectrographic 
measurements is limited to about | part in 10°, and thus is an error of about 
0-01 mass units (M.U.) or 10 Mev. in the region of heavy and medium- 
heavy nuclides. Since the radioactive decay energies are limited to a few 
Mev. only, the mass spectrographic data are not of much help in the under- 
standing of the radioactive decay of medium heavy and heavy nuclides. In 
recent years Duckworth etal. (1951) has reported some new atomic mass 
measurements in which the accuracy in the region of mass numbers 120 to 200 
is claimed to be of the order of 0-001 M.U. or about | Mev. 


Several attempts have been made since 1935 to set up a suitable 
semi-empirical mass-formula for the estimation of atomic masses when 
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experimental data are lacking and they have proved of immense help in the 
understanding of many nuclear properties. 


In the last few years more than 100 alpha-emitters have been discovered 
in the region extending from bismuth (Z= 83) to transuranium elements 
upto Z = 98 and in the rare-earth region (Z = 59 to Z= 72). The systematics 
of alpha decay properties have been well defined in a paper by I. Perlman, 


A. Ghiorso and G. T. Seaborg (1950) and the properties of possible alpha- 
emitters have been predicted. 


A number of short-lived alpha-emitters have been observed in the iso- 
topes of gold (Z= 79) and mercury (Z= 80) (Thompson e7 al., 1949). Of 
great interest is the discovery of a number of alpha-emitters in the rare-earth 
region with half-lives ranging from a few minutes to a few days and alpha- 
particle energies in the range 4-2 to 2:0 Mev. (Thompson et a/., 1949), Weaver 
(1950), Rasmussen et al. (1950 & 1953). 


A good test of the reliability of the semi-empirical mass-formula is to 
calculate the alpha-disintegration energies of a number of nuclides and com- 
pare them with the experimentally observed data. Several attempts have been 
made in this direction and the inadequacies of some of the semi-empirical 
formule in the investigation of the radioactive decay energies have been well 
established but the position with regard to the correction-term to be added to 
the semi-empirical formula so as to give better fit with experimental results 
in the regions of medium-heavy and heavy nuclei has not been clarified. 


In what follows, it has been shown that Fermi’s semi-empirical formula, 
with a suitable correction-term which depends on both A and Z, gives reliable 
estimates of alpha-disintegration energy in the medium-heavy and heavy- 
nuclide regions except for the magic number nuclei containing 126 neutrons 
but even in this case the agreement is better. The correction-term has been 
guessed with the help of Duckworth’s new atomic mass-values. 


Calculation of alpha-disintegration energy (E,) 


Let us consider an alpha-emitter of atomic mass M (Z, A) and the 


daughter atom of mass M (Z — 2, A — 4). The alpha-disintegration energy 
E, is given by 


E,=M(Z, A) — M(Z — 2, A — 4) — M(2,4)....M.U. (1) 


According to Fermi’s semi-empirical mass formula (Nuclear Physics, 1951) 
M (Z, A)= 0-99391 A — 0-00085 Z + 0-014 A2'3 + 0-083 x 


ID __ 2 
(A/2 — + 0-000627 Z?/A"? + 8 (A, Z) 
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where 
5(A, Z)=0 for A odd, Z anything 


— + (0-036 A-®4 for A even, Z lod (2) 


substituting (2) in (1) E, in Mev. is given by 


E, = — 27:90728 + 13-03610 sA?'® — (A — 4)?'3} — 77-28545 x 


(A — 22) + 0-58383 { © ». &-2 


ra oe ks + 031: 
A (A — 4) Alls ' A ays} + 93] 15 x 

x {8(A,Z) — 8(A— 4,Z — 2)} (3) 
The last term in (3) is zero for A odd, Z anything and = + 33-52140 


ae a _— 4\-3/4 _. A-~314 (even 
x (A — 4) A-3'4) for A even, Z lodd 
(1) Heavy nuclides.—The formula has been used for the calculation of 
E, in the isotopes of the heavy elements for Z = 96 to Z = 83 and the results 
are given in the Table I. 


It is seen that in every case the calculated E, comes lower than the 
observed one. For isotopes of a given element, E, is found to increase 
linearly with the decrease in the mass number A, the rate of increase being 
nearly the same for isotopes of different elements, i.e., it is independent of Z. 


For a given A, the calculated E, is found to increase with the increase 


of Z and*the rate of increase is nearly the same for all values of A ranging 
from 197 to 246. 


The observed E, shows a general trend of increasing with decrease of A 
(except in the neighbourhood of completed neutron and proton sub-shells 
where there is a sudden drop in E, on the lower side of A containing less 
than 126 neutrons [Perlman ef a/. (1950)]). The rate of increase depends on 
A and gradually decreases for higher values of A and Z. 


The difference E,ons. — Eacai. has been plotted against A. It is clear 
from the graph that this difference at first decreases slowly with increas- 
ing A but for A & 209 it shoots up abruptly—a maximum occurring for 
A=211, 212 and 213 for Bi, Po and At respectively (region for N= 126 
—completed neutron-subshell). Then it decreases rapidly as A increases. 
In the transuranium elements, the difference is of the order of 1 Mev. or 
even less. 


(2) Medium-heavy nuclides—There have been observed a number of 
short-lived alpha-emitters in highly neutron deficient isotopes produced on 
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TABLE [| 





ee, 


Observed Estimated 


aes Eq with 
slides a -e Estimated ; a 
Nuclides _—_ Half-lives EaXy 4 i in Mev. Difference the correc- 


in Mev. tion-term 
in Mey, 


cm . ' 
Cm**4 2 “00 
Cm*43 a: Pa 
Cm2!2 150d -28 
Cm* 55d ~ 
Cm 4° 26-8d -48 
Cm238 2-5h 61 


5-894 
196 
*336 
»* 465 
“619 
990 
998 


5-880 
030 
-173 
312 
-460 


+372 
534 
699 
875 
-062 
346 
633 


-449 
-738 
“118 
-307 


729 
-020 
-163 
*319 
461 

606 
-750 
927 
-044 
193 


130 
285 
-440 
648 
6-798 
6-961 


AAANAAnANH & 


gsAm 742 400) 
Am?!! 590» 
Am?!0 50h; 53h 
Am239 | 2h 
Am238 1-5h 


ora™ ~10¥ 
Pus ~6000)" 
Pure 24000) 
Pu238 92) 
Pu2s6 2 ‘Ty 
ee 8°Sh 
ra 22m 


t 
on ‘ 
= 

t t 


ogsNp??? 2-2 x 10%y 
Np?*5 435d 
Np238 35m 

53m 


oo” -51 x10%y 

2-46 x 107 

-91 x 108) 

+35 x 105) 

1-6 x 10°y 

70y 
4-2d 
20-8d 
58m 
9-3m 


Darn & HDRDBMAaan 
ARR ANFHHADH UUNNSSE 


ANUnLpHLPHF 


3-43 x 104) 
17-7d 
Pa**® 1-5d 
Pa228 22h 
ara 38m 
Pa226 1-7m 


DAM QaVnnnaaannk DAWN ADANMMAAN 


HHAHS EH UDHHHAHRHAWY 
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Nuclides Half-lives 


1-39 
8-0 104) 
7000) 
1-90)» 
18-6d 
30:9m 
7-8m 
Very short 


6: I3y 
21-Ty 
10:0d 

2:9h 

2:2m 

Short 


1622y 
3°64d 
11-2d 
38s 
31s 
Short 


4-8m 
27°5s 
‘025 
very short 
19m 


3°83d 
54°58 
3-925 
019s 


Em2! 9 
Em2?!8 
Em2!? 10-35 
Em246 ~—- Very short 
Em?!? 23m 


ggAt*t® 
At21? 
At?16 
At?15 
At?!4 
At2is 


‘021s 

l 0-35 

10-4 
Very short 


101°), 


Several secs. 


TABLE I (Continued) 


Observed 
A 
E. XA —4 
in Mev. 


Estimated 
E, in Mev. 


PRRRWHWWW 


Ap RWW WwW Ww 


WWW Ww WwW 


BWWWNNN PWCWWwW 
t 


WuwnNN bd 
t 


Difference 


NNNN = — wmwNNN— 


WwWwwNhva- 


NA PR HWW beet se oe 8 tap 


Dannie sh 


Estimated 
Eq With 
the correc- 
tion-term 
in Mev. 


*427 
-720 
*853 
-048 
5-203 
"3350 
506 
658 


-425 
*579 
*889 
043 
-192 
343 


*128 
-430 
‘576 
*727 
5-878 
030 


*254 
-410 
562 
*720 
650 


401 
*792 
-992 
*198 
394 
599 
395 


-530 
684 
816 
968 
-120 
*272 
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TABLE I (Continued) 


Observed Estimated 


—_ E, with 
Nuclides  Half-lives E,> —_ gry Difference _ the correc- 
! Re 


in Mev. tion-term 
in Mev. 


7°52 
6-12 
unobserved 
5-76 
STP 
5-85 
unobserved 
6-02 
6:32 
6°13 
6-48 
6-68 


wa 


-424 
- 568 
-719 
-879 
-022 
-176 
-322 
-477 
“648 
-774 
-960 
-169 


AFPHHLHAHHAHWWWWY 


NYADADDAAAWMUWN 


-23 
-02 
-64 
97 


-798 
-131 
*299 
-461 
- 589 
‘778 
-922 
-099 
*243 
346 
-495 
°732 
*882 
041 
-218 
341 
-483 
-648 


12 
57 

5-5] 
12 
34 
-30 
-40 

5-32 
56 
-67 


ADAMAAnAAnAUnAanpaRAPH HW 


ON 


82 
-96 


Se SSS KH BE QUNES 


Pesee yey ey yenyyyl > 
wr 


19-7m -72 
47m -09 
60:5m -32 
2-16m _ 6°87 
5-0d -95 
Very short “21 
62m 61 
25m -69 
9m -95 
2m *43 


802 
3-964 
-127 
-293 
-456 
‘647 
-028 
-178 
6-328 
6-479 
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211 216 221 22 23 
Mass Number A =" 


GRAPH 


irradiation of gold with high energy deuterons and these have been tenta- 
tively assigned to be isotopes of gold and mercury having mass numbers 
between 185 and 188 (Thompson and Ghiorso, 1949); Rasmussen, 1950; 
Rasmussen, Thompson and Ghiorso, 1953). The observed alpha-energy 
in case of gold is 2-20 Mev. and 5-7 Mev. for mercury. These values have 
slightly been revised but the masses of alpha-emitters have not been deter- 


TABLE II 


Observed Estimated 


A Estimated Eg with 
Nuclides _ Half-life EaXq_4 E,in Mey. Difference thecorrec- 
sa Rilo tion-term 

: in Mev. 


7gAul® wi = 4-058 . 5-888 
Aut8é : 18% -960 ‘ 5-795 
Aus? ron wie 3-868 ]- 5-708 
Autss - - 3-776 5-621 


solig ie -592 . 6-476 
He 7 -453 6-333 


a . 6-182 
— = en 





* Masses have been tentatively assigned. 
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mined conclusively (Rasmussen, Thompson and Ghiorso, 1953). Table [I 
given above shows the observed and calculated E, and their difference, 


The calculated E, is in every case lower than the observed one and the 
difference ranges from 1:2 to 1-7 Mev. which is nearly one-third of the 
calculated values. 


Inadequacy.—The genera! conclusion is that Fermi’s formula, though 
giving reliable estimates in the transuranium elements, is inadequate in other 
regions. The inadequacy of Bohr-Wheeler (1939) formula and of Bethe- 
Weizsicker (1936) has already been established by other workers (Saha and 
Saha, 1946; Das, 1950; Pryce, 1950). The empirical correction-term due 
to Stern (1949) which is applicable to A> 208, brings about better numerical 
agreement for heavy nuclides (Jha and Dube, 1952) but it is of no help in 
all the cases Z= 83, 84, 85 in the medium heavy and heavy regions. 


The failure of Fermi’s formula may perhaps be attributed to the absence 
of a suitable correction-term similar to the one proposed by Stern for heavy 
elements. 


Correction-term to Fermi’s formula 


Duckworth’s recent mass spectrographic data (Duckworth, Preston, 
1951; Duckworth, Strasford et a/., 1951) and their comparison with calcu- 
lated masses from Fermi’s formula [equation (2)] have made it possible 
to estimate the correction-term which is found to fit bettter with the observed 
values. The calculated masses from Fermi’s formula have been compared 
with the observed data for Z = 92, 90, 82 and 78 in the table III. It is seen 





TABLE III 








eee Mass calculated 
; Mass observed Mass calculate : : 
Nuclides sean oo ated Difference using correc- 
tion-term 
Pm a 238-1241-410 238-1205 0036 238-1218 
U2 235-1156-+-10 235-1130 -0026 235-1165 
U2s 234-1129-+-10 234-1107 -0022 234-1139 
eo 1 n?** 232-1093-+10 232-1069 -0024 232-1085 
aero *** 208 -0422 4-15 208 -0591 — +0167 208 -0436 
Pb2” 207 -0412 +15 207-0576 —-0164 207 -0439 
Pb20¢ 206 -0394 +15 206 -0550 -0156 206 -0425 
ght! 96 196-0274 +6 196-0379 -O105 196-0244 
pt % 195-0265 195-0366 —-O101 195-0245 


res 194-0256 -+-14 194-0341 — -0085 194-0232 
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that in some cases the calculated mass is higher than the observed one and 
in some cases lower, the difference 4 m ranging from + 0-004 to — 0-017 
mass units in the region 194< A< 238. 


From the isotopes of U, Th, Pb, Pt, it is clear that 4 m is a function 
of A and decreases with the decrease of A. But the dependence of correction- 
term on A only cannot be regarded as satisfactory. It depends on Z and 
increases with the decrease of Z. (It may be mentioned here that Stern’s 
correction-term is independent of Z and increases with A like this correction- 
term.) 


In the region of heavy nuclides, the experimental data are meagre. The 
atomic masses of Pt isotopes and of Pb2°* are known accurately. The masses 
of Pb2°? and P2°6 are determined from the known mass of Pb? and the 
observed neutron binding energy in Pb? and in Pb.?0 


M (207, 82) = M (208, 82) — Mn + Bn (208, 82) (4) 
M (206, 82) = M (207, 82) — Mn + Bn (207, 82) (4) 
where Mn= mass of the neutron= 1-00898 M.U. and the observed neutron 
binding energies in Pb? and Pb?®, i.e., Bn (208, 82) and Bn (207, 82) are 


7:38 Mev. and 6-73 Mev. respectively (Harvey et al., 1951, Kinsey et al., 
1950 and 1951). 


It is difficult to get a simple correction-term which may satisfactorily 
represent the correction for all ranges of A. In the region of medium heavy 
and heavy nuclides correction 


4 m=0-0290 — 3 [KA!:? — K’Z}°] in M.U. (5) 
where 
K=0-0002 and K’=0-00014. 
The masses calculated with this correction-term differ from the observed 
ones by values ranging from -0003 to -003 mass units in the medium heavy 


region. 
Calculation of E, using the correction-term 


With the introduction of this correction-term in Fermi’s formula, the 
calculated E, is modified considerably. The correction to E, in Mev. is 
given by 
4 E,=K, [Z'* — (Z — 2)'5]— K, [At:? — (A — 4)!-?] Mev. (6) 
where 

Where K,=0-39108 and K,=0-55869 for medium heavy and heavy 
nuclides. The calculated alpha-disintegration energy using this correction- 
term, i.e., E, + 4 E, are given in the last columns of Tables I and II. 
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DISCUSSION OF RESULTS 





There is a numerical agreement between the observed data and the 
calculated values using the correction-term in the region of heavy nuclides, 
In the region of 126 neutrons, there is some improvement, the difference 
between the observed and calculated values for At*!*, Po*!* and Bi*" is of 
the order of 4 Mev., 4-4 Mev. and 2-2 Mev. respectively which can satis. 
factorily be explained on the nuclear shell-structure. 


The mass numbers of gold and mercury isotopes which are short-lived 
alpha-emitters have been tentatively assigned to lie between 185 and 
188 and the alpha decay energies are 5-07 and 5-6 Mev. respectively. 
These results agree with the calculated ones for the isotopes having mass 
numbers 188. 
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KINETICS OF THE OLEFINE-BROMINE REACTION 


Part VII. The Reaction between Cortonic and Dimethyl Acrylic Acids and 
Bromine in Carbon Tetrachloride* 


By I. M. MATHAI AND S. V. ANANTAKRISHNAN, F.A.Sc. 
(Department of Chemistry, Madras Christian College, Tambaram, S. India) 
(Originally received January 19, 1954 
Received in revised form May 27, 1954) 


Srupies of the reaction in glacial acetic acid have been reported in earlier 
parts.* It was then considered desirable to extend the observations to non- 
hydroxylic solvents and preferably nonpolar solvents. The results of experi- 
ments in carbon tetrachloride are reported here. 


From the trend of observations in other solvents» * that had been 
reported by different workers and from observations of the reaction in the 
gas phase,® it was anticipated that the reactions in this solvent will involve 
a polar surface. Now, all heterogeneous reactions involve the adsorption 
of one or more reactants in a single or a multimolecular layer and it is neces- 
sary also to note the influence of the nature of the surface on such adsorp- 
tion. In the exploratory work, it is desirable to avoid other complications 
and in spite of the disadvantages of choosing carboxylic acids, these com- 
pounds have been taken up for study. This avoids the complication of the 
peroxide effect known to be negligible in these and also enables a comparison 
with the work in other solvents. 


EXPERIMENTAL 
(i) Material: used 





Carbon tetrachloride-—Pure carbon tetrachloride was kept in contact 
with an acid solution of potassium dichromate for three days with frequent 
shaking. The carbon tetrachloride was then separated, washed with distilled 
water, then with a dilute solution of potassium hydroxide and then again 
with distilled water. The moist solvent was then dried over anhydrous cal- 
cium chloride for several days and fractionated. The middle fraction, 
b.p. 76:6-76-8/760 mm., was collected and stored over sticks of potassium 
hydroxide in amber coloured bottles. The use of phosphorous pentoxide 
for drying this solvent was found to be unsatisfactory and gave erratic results 
for the course of this reaction. A report of the dielectric constant studies 
on this apsect of the work is being published separately. 


* Based on the thesis of ILM.M. accepted for the M.Sc. degree of Madras University. 
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Trans-Crotonic acid.—The product supplied by Tennessee Eastman 
Corporation was used after purification by repeated fractional crystallisation 
from petroleum ether (60-80°) to constant m.p. and stored in the dark, over 
anhydrous calcium chloride in a vacuum desiccator. The choice of this 
procedure was necessitated by the need for minimising the chances of con- 
version to the cis form. Freshly recrystallised samples, m.p. 72° C., was 
used for each experiment. 


2:2 Dimethyl acrylic acid—The compound was prepared essentially 


by the method of Barbier and Glaser.’ Mesityl oxide (1 mole) was oxidised 
by an alkaline hypochlorite solution (2 moles hypochlorite and 4 moles sodium 





80 80 


Per-cent. Reaction 











600 
Time in minutes 
Fic. 1 


Reaction between dimethyl acrylic acid and bromine at different concentrations. 
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hydroxide) keeping temperatures below — 5°C. Excess hypochlorite 
was destroyed by bisulphite and after the removal of the byproduct, chloro- 
form, the solution was acidified with dilute sulphuric acid, keeping tempera- 
tures low. Dimethyl acrylic acid separated on standing and a further yield 
was obtained by ether extraction of solution (yield 70% of theory). The 
acid was recrystallised from hot water and stored in a desiccator in the dark 
over phosphorous pentoxide, m.p. 70° C. 
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Bromine was prepared by the method indicated by one of us earlier and 
the purity checked from time to time.! 


(i!) Materials added for increasing surface 


Sodium bromide, barium sulphate and lithium bromide, of analytical 
reagent quality, were dried at 120°C. and stored in a desiccator before use. 
Alumina, silica and pyrex glass beads were used to provide additional surface. 

A4 
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The beads were treated by the same procedure, with reagents used at the 
same concentration and in the same order and prepared fresh for each experi- 
ment. The beads were soaked in a chromic acid solution for about 12 hours, 
with frequent shaking, washed thoroughly with distilled water, refluxed 
thrice with distilled water for 10 minutes at a time and then dried in an air 
oven at 120° C. for three hours. The beads used were uniform in each case, 
the pyrex beads being spherical and the others cylindrical. As a first approxi- 
mation, the surface area was calculated for these shapes from the dimen. 
sions and in the different runs, the surfaces were assumed to be proportional 
to the weights. The relevant data for the purpose are given in Table I. 


TABLE | 





Average radius 
” Average 





Material of bead ennes 1 cues | — 
External Internal sq. cm. | & 
: 
Pyrex oe] ve 0-287 cm. 1035 377 
Silica ..|  0+260 0-125 | 0-675 | 1-959 456 
Alumina oe 0-293 0-157 0-638 | 2-188 | 260 





(iii) Reaction Vessels and Measuring Vessels 


250 ml. and 100ml. amber coloured Pyrex bottles with well ground 
stoppers were used as the reaction vessels. The bottles were kept with chromic 
acid for 24 hours, washed thoroughly with distilled water and then steamed. 
The vessels were then dried in an air oven at 120°C. for about four hours 
and kept overnight in a vacuum desiccator. 


All measuring vessels were standardised and coated with black enamel 
and were submitted to the same treatment as the reaction vessels and were 
used only after 24 hours at room temperature. Under these conditions 
volume changes were found to be negligible. 


The reactions were carried out in an electrically maintained thermostat 
(temperature controlled to within 0-005°C.). All vessels were tested for 
evaporation losses, using bromine solutions. The course of the reaction 
was followed by the iodometric estimation of the unreacted bromine, in 
aliquot parts, from time to time. Rapid estimations were necessary as some 
of the dibromides formed slowly liberate iodine. 


To economise space, only a few typical runs are presented in the 
accompanying tables and figures. 
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TABLE | 
he Reaction between Bromine and Trans-Crotonic Acid at 35-1° C. 
ri- aed as =o 
rs, Concentration of | 0-2000 M 0-1000 M 0-0590 M | Q-0500 M | 0-0500 M 
d olefine | | 
a Concentration of 0-2093 M 0-1004 M | (+0508 M | 0-05055 M | 0-0503 M 
ir bromine A B | C | D E 
) | 
$e, - a 
xl A Time (min.) -- 0 2:25 7-25 14-5 25-5 45 83-25 350 [1086 2300 
Titre .. 19-60 19-60 19-50 19-20 18-80 18-20 17-80 16-25 | 13-90 11-50 
N- % Br reacted 0 0 0-51 2-04 4:08 7-15 9-18 17-09 | 29-08 41-32 
| 
al B Time (min.) -- 0 38-0 94-11 199-4 783 2705 5465 8102 | 
Titre -- 9-40 9-40 9-10 9-00 8-50 7-40 6-35 5-40 | 
% reaction -+ O 0 3°19 4-26 9-57 21-27 33-21 42-56 
C Time (min.) .. 0 8-05 215-5 744 1114 2749 5265 6768 = | 
— Titre -- 9-20 9-20 9-10 8-80, 8-55 7:70, 6:55 5-95 | 
% reaction ++ 0 0 1:09 4+35 | 7-07 16+30 | 29-47 34-25 
= D Time (min.) -- 0 106-8 227-1 683 1060 2494 3221 3983 Pyrex Bead 
, Titre -- 980 9-00 8-70 7-95 7-40 5-45 4-90 4-50 Surface 
% reaction .. 0 3-23 6-45 14-32 20-43 41-40 47-32 51-60 | 585-2 sq. cm, 
E Time(min.) -- © 136-25 312 1051 1780 2517-5 3832 5991 
Titre +» 9-20 9-00 8-85 7:85 | 7-00 6-40) 5-35 4-00] Silica bead 
% reaction ++ O 2-17 3-80 14-68 23-91 30-45 | 41-85 53-26 
nd TABLE II 
ric 


4 Reaction between Bromine and 2-2 Dimethyl Acrylic Acid at 35-1°C. 




















rs | | ’ 
Concentration of olefine | 0-1000 M 0-0500 M 0-0250 M 
| 
Concentration of bromine | 0-1009 M | 0-04956 M 0-02532 M 
1€ | | 
A | B 
re | | . | 
ns = es a ae : ay 
A Time (min.) sa 0 1-08 3°62 6°17 12-33 16-70 | 23-40 
‘litre (ml.) ++ 19-50 18-00 16-20 3°75 10-80 9-15 | 7-10 
% Bromine reacted .. 0 7-69 16-42 29-49 44-62 53-08 | 63-08 
‘at —_ ——_—- 
. Bb Time (min.) ce 0 11-70 21-43 29-07 36-53 42-74 61-50 
or Titre we] 9°35 9-35 9-15 8-40 7:40 6-45 4-90 
on % Bromine reacted .. i) 0 2-14 10-16 | 20-85 31-01 | 47-62 
; | — 
in C Time (min.) *- 0 328 687 775 | 850 942 985 (999 
ne Titre (ml. ) .- 10-00 9-70 9-30 8-85 | 8-20 6-46 4-95 | 4-00 
% Bromine reacted .. 0 3-0 7-0 11-5 18-0 36-0 50-5 60-0 
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TABLE IIT 


Influence of Surface on the Reaction between Bromine and Olefine for an 
Olefine Concentration of 0-0500 M at 30-1°C. 


Bromine Concentrations: same order but slightly different. 
A: Typical runs with added pyrex beads; B: Typical runs with silica beads. 


Crotonic acid Time ee 106-8 | 272 683 1060 | 2494 | 3983 
Titre ..| 9°30 | 9-00 | 8-70 | 7-95 | 7-40] 5-45 | 4-59 
% reaction wok 3°23 6-45 | 14-32 41-40 | 51-60 





2-2 Dymethyl Time oe °7 19-50 | 31-37 
acylic acid Titre ar +4 -0 9-10 | 8-40 
reaction . 3°85 | 12-50 | 19-23 


= 
y — 





Crotonic acid Time Pe | 1051 1780 
Titre eo| 920 | : 7°85 | 7-00 
% reaction = | Zs 14-68 | 23-91 | 36- 53-26 


wo 
~~ 





Per-cent. Reaction 


2-2 Dimethyl | Time J -22 | 23-87 | 29-97 | 39-52 | 46-95 | 7 
acry lic acid Titre NP -55 | 9-85 | 8-20 *25 | 6-40 | 
% reaction a . | 15- | 22-47 3i1- | 39-42 | { 
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DISCUSSION OF RESULTS 


Schwab® has remarked that where there is a surface reaction between 
ethylene and the halogens, both reactants are slightly adsorbed and the experi- 
ments of Gwyn Williams® have borne out these. The typical runs are reported 
below and the illustrative curves (Figs. 2, 3 and 4) show the trend. 


An inspection of these shows clearly the presence of an induction period 
which is sensitive to changes in the concentration of reactants as well as 
changes: in the structure of the olefine and in the case of crotonic acid the 
reaction even appears to have a negative temperature coefficient. It is gene- 
rally considered that reactions in carbon tetrachloride simulate a gas reac- 
tion. Preliminary experiments have shown that both reactants are adsorbed 
in varying amounts. A fuller study of this is in progress. The hetero- 
geneous reaction is unmistakably shown by the elimination of the induction 
period and a quickening of the reaction by packing the vessel with beads 
of different materials, 
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A qualitative comparison shows that the effectiveness of the surface 
in catalysing the reaction is in the order 































Pyrex glass > Silica (fused) > Alumina (fused) th 

A rough comparison can be made by transforming the results obtained so tl 
that they refer to the same surface area. As a first approximation, one may a 
assume that the time required for a given fraction of the reaction is inversely a 
proportional to the surface. The results obtained by such computation are re 
given in Tables IV and V below. if 
si 

TABLE IV 






































Reaction between Dimethyl Acrylic Acid and Bromine: at Constant 0 
Surface Area , 
oe ieee 
Time for 25% reaction Time for 50% reaction it 
Concentration Surface - oe Een - r 
sc | wc | wc | 35° C, P 
— s 
Pyrex oa 16-1 40-0 
0-066 M Silica wa 30-8 80-9 | 
t 
Alumina = 39-4 86-7 
= i. 
Pyrex ar 42-6 17°3 70°8 30-4 r 
0-05 M : 
Silica a 53-9 25-3 111-2 | 4855 I 
| Pyrex ‘is 139-0 216-0 ¢ 
0-025 M ¢ 
Silica on 197-0 283-0 
S 
Note: There have been small variations in the bromine concentrations. 
f 
TABLE V 
1 
Reaction between Trans-Crotonic Acid and Bromine at Constant Surface { 
Area at 35°C. = 
] 
Concentration Surface Time for 10°¢ reaction Time for 25% reaction Time for 50% reaction s 
min, min. min. A 
—_ a ote i i 
Pyrex “ 450 1103 4622 ( 
0-1 M 
Silica oe 665 3036 § 
Pyrex a 535 1172 4207 
0-05 M < 


1984 


Silica 
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An inspection of the tables clearly shows pronounced differences between 
the two acids. While consideration of structural factors may be deferred, 
the first stage in any mechanistic study is an analysis of the reaction order 
and the trend of rate constants. In spite of the choice of a nonpolar solvent, 
a heterogeneous reaction in solution may not show the simplicity of a gas 
reaction. The present observations are also inadequate for assessing the 
influence of the products of the reaction. It is also likely that there is a 
simultaneous homogeneous reaction. 


Taking the time required for a given fraction of a reaction, the figures 
obtained lead to an approximate relationship 


(C,/C2)? = ta/t 


in the initial stages, but, the right-hand side tends to smaller values as the 
reaction proceeds, the equality being nearly reached when the square of the 
concentration terms are used. With silica surfaces, the third power relation- 
ship does not hold even at the start and falls off later to the first power. 


In acetic acid solutions, Robertson and co-workers’ find that the reac- 
tion is of the third order at the low temperatures at which they worked while 
Anantakrishnan and Venkatarman working at or above 30° C. find that the 
reaction is of the second order. The gas reaction according to Gwyn Wil- 
liams® is apparenily of the third order. The difficulty in accepting a third 
or higher order for the reaction is the definite positive temperature coeffi- 
cient, though the position is rather complex with crotonic acid which has a 
small negative temperature coefficient. 


It is quite clear, however, that the reaction requires a polar environment 
for initiation of the reaction, though, it is just possible, that the polarised 
molecule after leaving the surface may initiate the homogeneous phase of 
the reaction in the nonpolar solvent also. The data so far obtained is not 
adequate for the purpose of separating the two phases of the reaction. 
Preliminary experiments have shown that the addition of dry crystals of 
sodium bromide or barium sulphate accelerate the reaction while the relative 
influences of the three different surfaces studied indicate that the presence 
of both cations and anions at the surface facilitate the reaction. Using only 
glassy surfaces, it is not practicable to study the geometrical factor but com- 
plex silicates appear to afford a better environment than either silica or 
alumina. 


We have to thank the Tennessee Eastman Corporation for a gift of the 
crotonic acid used. 
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SUMMARY 


The reaction between trans-crotonic acid as well as 2: 2-dimethyl acrylic 
acid and bromine has been studied in carbon tetrachloride solution in the 
dark. The reaction is found to be quite a complex one involving both homo- 
geneous and heterogeneous phases. The course of the reaction is governed 
by the nature of the olefine and the nature of the catalysing surface. 
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